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Reversible reaction-diffusion systems display anomalous dynamics characterized by a power-law 
relajcation toward stationarity. In this paper we study in the aging regime the nonequilibrium 
dynamical properties of some model systems with reversible reactions. Starting from the exact 
Langevin equations describing these models, we derive expressions for two-time correlation and 
^^ I autoresponse functions and obtain a simple aging behavior for these quantities. The autoresponse 

fvj . function is thereby found to depend on the specific nature of the chosen perturbation of the system. 
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I. INTRODUCTION 

The intriguing aging processes observed in nonequilibrium systems with slow (i.e. nonexponential) dynamics have 

r^ ■ been the focus of many intensive research efforts in the past. Whereas the initial studies almost exclusively focused 

O ] on systems like glasses [i| and spin glasses [3] , it has been realized quite recently that aging processes do not uniquely 

^5 . characterize these complex systems, but that they are also encountered in much sirnpler situations I3|. Thus, the study 
of aging phenomena has for example been extended toward magnetic systems (J, la, |6| , prepared initially in a disordered 

^ high temperature state and then quenched to or below their equilibrium critical point, and toward reaction-diffusion 

Ci • systems quenched to their nonequilibrium critical point |7| . These investigations of simple systems displaying aging 

c/j I have led to an increased theoretical understanding of aging processes taking place far from equilibrium. 

^ ' Diffusion-limited irreversible reactions are often characterized by the presence of an absorbing phase transition, 

2 ] separating an active phase from an inactive or absorbing state from which the system can not escape [3, S ll3] • These 

^ . noriequilibrium phase transitions have attracted much interest, and different universality classes have been identified 

' , ' [Sl, lll|. Various quantities, as for example the density of particles, display simple power-laws when approaching 

Zh . these nonequilibrium critical points. In addition, the dynamical correlation length increases as a power-law of time, 

Q ' similar to what is observed at an equilibrium critical point, revealing the presence of slow dynamics. Consequently, 

^ ' aging pro cesses have been studied in irreversible reaction-diffusion systems quenched to their nonequilibrium critical 

'""'■ point [i3, [ij, H) \lM, Mj [13, [13- Interestingly, two-time quantities, like the autocorrelation function C{t,s) and 

. ] the autoresponse function R(t,s), display in these systems the same simple scaling behavior as the corresponding 

^ . quantities in equilibrium critical systems: 

^; C{t,s) - s-'^fcit/s) (1) 

a^■ R{t,s) ^ s-^-''fR{t/s)0{t~s) (2) 

.• ' where a and b are nonequilibrium exponents, whereas the scaling functions fc and //? only depend on the ratio t/s, 
S^ : with fc{y) '^ y^-^"^/^ and fR{y) ~ y-^ril^ ^^j. y ^ \ Here z is the dynamical exponent, and Ac and A^ are called 
(-^ autocorrelation and autoresponse exponents. Finally, the step function d(t — s) in the expression for the autoresponse 

f— ^ , function ensures causality. The absence of detailed balance in irreversible reaction-diffusion systems reveals itself 
J* ' mainly by the fact that the relation a = b, trivially observed at equilibrium critical points, is no longer fulfilled (see 
. !^ I [3] for a recent review containing a thorough discussion of this point) . 

^^ . All these recent studies investigated aging at nonequilibrium phase transitions. This is of course due to the fact 

^ ■ that in reaction-diffusion systems with irreversible reactions nonexponential relaxation is exclusively encountered at 
absorbing phase transitions. 

It is very remarkable that non-exponential relaxation is generically found in reversible reactions, and this without 
fine-tuning of the system parameters (as would be needed in order to be exactly at a phase transition point) . A power 
law behavior in the long time limit was first predicted in [19[ for the bimolecular reversible reaction A + B ^ C 
taking place in solutions, based on physical arguments involving spatial concentration fluctuations. This power law 
approach to stationarity in reversible reaction-diffusion systems was later verified through more elaborated approaches 
[20, 21, 22, 23], through numerical simulations 24, 25(], and through some exactly solved models [26, 27, 28, 29,]. 
In addition, this power law behavior has been observed experimentally in excited-state proton transfer reactions 
[30, [Sll, IsJ- ^11 these studies demonstrate that the most important ingredient for aging, namely slow dynamics, is 
typically encountered in reversible reaction-diffusion systems. 

We propose to extend the study of aging phenomena to these reversible reaction-diffusion models. Starting from 
idealized reaction-diffusion models, we exploit the fact that a set of Langevin equations, describing the time evolution 



of the system, can be derived exactly for these systems [26| . Within the standard field theoretical representation of 
reaction-diffusion models we derive exact expressions for the two-time correlation and autoresponse functions in the 
aging regime. For all the studied models we recover a simple aging behavior where the scaling functions of these 
two-time quantities only depend on the ratio of the two times. The most remarkable result obtained in this study 
concerns the autoresponse function whose expression is found to depend on the nature of the perturbation applied to 
the system. 

The paper is organized as follows. In Section |TT1 we consider the reaction A + A ^=^ C. Starting from the master 
equation description of the model, we use the exact map to a set of two Langevin equations for some random complex 
variables a and c. The concentrations na{t) and ndt) of the particles of types A and C are then given by the 
average of the random variables a and c over the complex noise: Ua = (a) and ric = (c). Exploiting the controlled 
approximation scheme established in [26[, we derive exact expressions for two-time correlation and autoresponse 
functions in the dynamical scaling or aging regime. In Section [llll we extend our study to the reactions A + B ^^ C 
and A + B ^ C + D. Using the same approach as in Section |TT1 we find the functional dependence of the two- 
time quantities also for these reactions. In Section IIVI we summarize our findings and discuss open problems. Two 
additional technical points are discussed in the Appendices. 

II. THE A + A^C REACTION SCHEME 

A. Model and exact Langevin equations 

Following [2y|, we consider two types of particles (called A and C particles) that diffuse on a d-dimensional hy- 
percubic lattice. Allowing multiple occupancy of a lattice site, particles at the same site may undergo the following 

reactions: (1) the reaction A + A -^ C with rate Aq where two A particles recombine to form a C particle and (2) 
the reverse reaction C ^ A + A where a C particle dissociates with rate /i. These microscopic rules are readily 
translated into the following master equation for the probability P {{m}, {n}; t) (where {m} ~ {■ ■ ■ , rn^, • ■ • } resp. 
{n} — {■ ■ ■ ; rij, • • ■ } are the occupation numbers of particles A resp. C for every lattice site i) that the configuration 
{m}, {n} is found at time t: 

|P ({m}, {n}; i) = ^ H.P ({m}, {n}; *) + E E ^^^^ (i™}' i"): *) • (3) 

« « Jii) 

Here the term TiiP ({m-}, {n}; t) contains the reactions taking place at lattice site i and is given by 

TiiP{mi,ni;t) = ^i[{ni + l)P{mi - 2,ni + l;t) - niP{mi,ni;t)] 

+Ao [{mi + 2)(m, -t- l)P{mi -)- 2, n^ - 1; t) - m,{mi - l)P(m,, n,;t)] (4) 

where for convenience we wrote as arguments of P only the occupation numbers mi and Ui of lattice site i that are 
changed by these reactions. Diffusion processes are captured by the second term in Eq. ^ where the sum over j{i) 
is a sum over the nearest neighbor sites j of the lattice site i. Indicating again only the occupation numbers that are 
changed in the process, we have the following expression for DijP: 

T^i.jP = -T^[{mj + l)P{mi-\,mj + l\t) -miP{mi,mf,t)] 

+ -^ [{fij + l)P{n, ^ l,nj + l;t)- niP{n„nj;t)] (5) 

where Da and Dc are the diffusion constants of the A and C particles and i is the lattice constant. 

The master equation ([3]) must be supplemented by initial conditions. We here consider the case of an uncorrelated 
Poisson distribution on each site and for each particle species. Prepared in this initial state, the system evolves toward 
chemical equilibrium in the long time limit |26j . 

As shown by Rey and Cardy |26| the dynamics of this model allows an exact description in terms of a set of coupled 
stochastic Langevin equations. Introducing A ~ Aq^'^, one gets 

{dt - Day^)a{x,t) = -2Aa2(x,i)-h2/xc(x,t)-f C(x,i) (6) 

(5t-DcV2)c(x,i) = Aa2(x,i)-/xc(x,t), (7) 

where C is a complex Gaussian noise with zero mean value whose correlation is given by 

(C(x, t)C(x', t')> = 2(^c(x, t) - Aa2(x, t))d{x - x')d{t - t'). (8) 




Here, the vector x describes the d-dimensional space coordinates, whereas the bracket notation stands for the average 
over the noise. Inserting Eq. ([7]) into this expression yields 

(C(x, i)C(x', i')) = -2dt{c{t))S{^ - ^')S{t - t') (9) 

for spatial homogeneous initial conditions. 

Note that the variables a{yi,t) and c(x, t) do not represent the particle densities, as they are complex |26i]. The 
mean densities na{t) and ndt) of the particles of types A and C, which are of course real valued, are given by the 
averages (a(x, t)) and (c(x, t)) of these complex variables over the noise. 

As the density nc{t) = (c(x, i)) reaches a stationary value in the limit t — > cx), the noise correlation, see Eq. 
([9]), will vanish in the long time limit. As a consequence of the vanishing of the fluctuations at equilibrium one can 
compute the actual values of the equilibrium densities, which are given by their mean field values. These equilibrium 
densities Ooo and Cqo satisfy the relationship 

Aa^ = MCoo , (10) 

which follows directly from Eq. ([8]). 

It is important to notice that the dynamics conserves the quantity na{t) + 2nc{t) due to the total mass conservation. 
In the case of equal diffusion constants Da — Dc = D, which is the case we discuss in the following, we have in addition 
that X == ci + 2c obeys the noisy diffusion equation 

(9t-I?V^)x(x,t)=C(x,i), (11) 

as follows directly from Eqs. ([6]) and ^. Since {() — and the initial conditions are homogeneous, the average value 
of the field (x(x, i)) is conserved. In particular we have xo = (x(x, 0)} — (x(x, oo)) = Xoo- Using Eq. (fTO|) one can 
therefore express the equilibrium densities as a function of xo [26,] : 

and Coo = -(xo -floo) ■ (12) 

Starting from this Langevin description and exploiting the existence of the conserved quantity, Rey and Cardy [26[ 
developed a systematic approximation scheme that enabled them not only to derive the power-law relaxation toward 
equilibrium but also to compute the corresponding amplitude exactly. In the following we extend this approach to 
two-time quantities and derive in leading order exact expressions for correlation and response functions in the aging 
regime. 

B. Two-time correlation function 

New insights into the behavior far from stationarity can be obtained through the analysis of the connected correlation 
function 

Cit, x; s, y) = (<5c(x, t) Sc{y, s)) - (5c(x, t)) (<5c(y, s)) , (13) 

with (5c(x, t) = c(x, t) — (c(x, oo)) = c(x, t) — Coo- In the definition of the correlation we have taken into account that 
in the stationary state the mean value of the variable c is Coo- In case x = y, we are dealing with the autocorrelation 
function C{t, s) — C{t, x; s, x) — C{t, 0; s, 0), where in the last identity we exploited the spatial homogeneity of our 
system. 

In principle, we should define a similar quantity for the variable a(x, i). However, as 5a(x, t) — a{'x.,t) — Ooo — 
(5x(x, t) — 2 ^c(x, t) with (5x(x, i) = x(x, i) — Xoo, we can immediately derive the correlator for a once we know the 
correlators for the conserved quantity x and for the variable c. 

Starting point for the computation of the correlation function is the following Langevin equation for the quantity 
<5c(x,i) [261]: 

{dt-DV'^ + a)Sc{x,t) = 4A5c2(x,t)-4A5x(x,t)(Jc(x,i)-fA(5x^(x,f) 

+ ^(^-M)'5x(x,t), (14) 

where a = -iXa^ + u. 



The formal solution of this equation is 

5c = G[{c - ca)5{t) + 4A Sc^ -AXSxSc + X 5x^ + (o" - A^) ^x] , (15) 

with the Green function 

G[f] (x, t) = / At' I A^x' e-^(*-*')Go(x - x', t - t')/(x', t') . (16) 



The nonlinear equation (|15p is solved by the systematic approximation scheme developed in [26|. As a result we 
obtain to leading order the following expressions for ((5c(x, i)) and ((5c(x, i)(5c(y, s)): 

(<5c(x,t)) = ^(5x'(x,i)), (17) 

{5c{^,t)5c{y,s)) = {^^^ {Sx{^,t)Sx{y,s)). (18) 

It follows that the calculation of our correlation function reduces to the calculation of the expectation values {Sx'^ (x, t)) 
and {Sx{^, t) (^xIY) *))■ R-oy and Cardy already computed the first quantity and obtained in leading order 

(5x'(x, t)) = 2(co - c„o) iSirDt)-''/^ . (19) 

For the computation of ((5x(x, i) i5x(y, s)) we exploit the fact that x(x, i) obeys the noisy diffusion equation pi[) with 
the noise correlator given by Eq. Q. This readily yields the expression 

(^X(x, t) <5x(y, s)) = -2 y clti y d''xiGo(x - xi, i - ti)Go(y - xi, s - ii)at (c(ti)) , (20) 

where Go(x, t) is the free propagator: 



Go(x, t) = eit) {AirDt)-'"^ exp (-^) , 



(21) 



and 6{t) is the step function. Hence 

(^X(x, t) <5x(y, s)) = -2 y' At, [47rD(< + s - 2t,)]-"^ cxp |-^^^|i-^L| a*(c(ii)) . (22) 

For X = y we obtain up to some numerical prefactor an integral that Rey and Cardy already discussed in [26| . In the 
aging regime, where both t and s are large, this then leads to the following leading behavior: 

(<5x(x, i)^x(x, s)) - 2(co - c^) [^nD{t + s)]-"'^ . (23) 

We therefore have for the autocorrelation function in the aging regime the expression 
G(i, s) = ((5c(x, t) (5c(x, s)) - ((5c(x, t)) ((5c(x, s)) 

'''^^ {6x{^,t)Sx{y,s))-(^] (<5x'(x,t))(<5x'(x,5)) 



2cr 

= K^) i^o-Coo)[^7TDit + s)]-''/' , (24) 

where in the last line we omitted all sub-leading correction terms. It is interesting to notice that the term 
((5c(x, t)) {Sc{x,s)} is of the order g-d/2f~d/2 ^^^ therefore only contributes to the correction terms. This reflects 
the fact (already noticed in [20|) that the quantity c is not a self-averaging quantity. 

The derived expression for the autocorrelation function can be cast in the form C(t,s) ~ s'^^fc{t/s) that char- 
acterizes a simple aging behavior. This also yields the following values for the nonequilibrium exponents (see Eq. 

©): 

b = d/2 and Xc/z = d/2 . (25) 

The space-time correlation function C{t,ii.] s,y) can be computed in exactly the same way, with the final result 

C(i,x;s,y) = i(^^) [c,-c^)[A^D{t + s)]-''/^e^j>^-^0-^ . (26) 



C. Two-time response functions 

The system can be perturbed in several possible ways in order to compute the linear response. One of the possibilities 
is to inject new particles (which can be particles of either type A or C) at time s. The 'injection' process is assumed 
to be random with the same small occurrence probability at each lattice site. The response of the system to that 
perturbation is then monitored at a later time t by measuring the densities of particles of type A or C. In this way we 
obtain different responses that we note as R( {t, s) where i stands for the type of particles that are created whereas / 
indicates the type of particles whose density is measured. Thus, the response R^it, s), formally given by the equation 



,Au „X _ ^(«)W 



R^it,s) 



Shc{s) 



(27) 

he — ►O 



means that we are measuring the linear response of the A particles density to the additional creation of C particles 
only. It is important to notice that this process violates the conservation of the quantity (x) ■ In order to assess the 
impact of this violation on the response of the system, we also consider a process that conserves the total mass of the 
particles. As we discuss in the following, the response of a reversible diffusion-reaction system strongly depends on 
the chosen perturbation. 

Let us start by injecting particles of type C into the system and by monitoring the subsequent change in the particle 

density of the same particle type. A creation process -S C modifies the single site part of the master equation ([3]) 
which now reads as 

7^,f P(m,, n,; t) = Ti^Pim,, nf,t) + hc{t) [P{m,,n, - 1; t) - P{m,,n,- 1)] , (28) 

where HiP{mi,ni;t) is the expression ^ one has without additional creation processes. 

We might want to inject Qc additional C particles at time s, so that the particle injection probability hc(t) is 
given by 

hc{t)^ncS{t-s) . (29) 

However, as the final state of the evolution will again be a homogeneous state, the concrete form of he is not of real 
importance as long as the particle injection has ended before the measurement of the response at time t, with the 
total number of injected particles being given by 

nc= [ drhcir) . (30) 

Jo 

With this particle injection process the Langevin equation for the C particles is now given by 

{dt - DV^)c''{x, t) = Aa'Xx, tf - /xc''(x, t) + h{t) , (31) 

where in the continuum limit we set 

Mt) = ^ and n = ^, (32) 

thereby dropping the index C for the quantities divided by the volume of the system. We use in the following the 
index h in order to emphasize the presence of the additional creation process and to distinguish the corresponding 
quantities from those obtained in the absence of this process. The Langevin equation for the A particles is unchanged 
by the creation of C particles, but the noise-noise correlator ^ becomes 

(C(x,t)C(x',t'))=2[Mt)-9,(c^t))]5(x-x')<5(t-t') , (33) 

for spatial homogeneous initial conditions, due to the Eqs. ^ and (PT|) . 

With this one readily verifies the validity of the following set of algebraic equations: 

at + 2c^ = ao + 2co + 2^! - xo + 2r! , (34) 

A(a^)' = A^4. (35) 

Note that as il has the meaning of the average total change in the number of C particles due to the creation process, 
Eq. (|34p reflects the modification of the total mass due to this process. Equations p4|) and pS]) immediately yield 
the following expressions for the mean values a^ and c^ in the new stationary state: 




aSo = ^K/l + ^(Xo + 2r!)-lj and cj, = l(xo + 2f7 - a^o) (36) 



which differ from the expressions p2|) without the field h through the replacement of xo by xo + 2^1. 

As the coefficients a^ and c'^ are again constants, it follows that also in this case the perturbation series derived 
by Rey and Cardy [2^ and sketched in the previous Section only involves time independent coefficients. Therefore, 
any quantity which is a function of the coefficients a^ and c'^ can be derived in exactly the same manner as discussed 
before. For instance, the constant a becomes 



4Aa^ 



A^ 



(37) 



Hence, in order to calculate the /i-dependent density (c'') = cj^ + {Sc^) we just need to substitute into Eq. P7)l the 
expression (j37]) and the /i-dependent expression for {5x^)- This last quantity is obtained by noicing that due to the 
presence of the additional injection process the Langevin equation for x ^^o^ reads 



^ + V2')x''(x,t)-C(x,i) + 2Mi) 



(38) 



This yields for {{5x Y) the expression 



h\2\ 



iiSx") 



2(co - c'^)i8nDt)-'^/^ + 2 /" dr [87:D{t - r)] 

Jo 



-d/2 



h{r) 



(39) 



in leading order. 

The last remaining step is to calculate the two-time linear response function 



R^it,s) 



S{cHt)) 



Sh{s) 



which after some elementary algebra is given by the expression 



Rcit.s) 



Xfj. 



2 r 



24:Xfi 



[Co 



+2 



^[8nD{t-s)]-'/'e{t-s) 



(40) 



inDty'^/^eit - s) 



(41) 



where we used the fact that 3-^777 = 0{t — s). This expression can also be cast in the standard scaling form R^it, s) = 
s-^^°-fjj(t/s)d{t — s), with a — d/2 — 1 and Xr/z — d/2. Interestingly, we have a j^ b, even though the system evolves 
toward chemical equilibrium. We also note that Xr = Xc- 

The change in the density of particles A due to the injection of C particles, i?^, is related to the response R^ by 

i?^(i,s)=2^(i-s)-2i?g(i,s) , (42) 

which follows directly from (a'*} + 2(c''} = (x'') = Xo + 2ri. Thus the expression (|42p is the sum of two terms with 
different scaling behaviors where the constant term is the leading one. 

Let us now discuss the response of the system to an injection of particles of type A. One possible way of doing this 

consists in injecting pairs of A particles into the system with a rate h2A'- — > 2 A. As these pairs can immediately 
react to form a C particle, it is expected that this leads to the same behavior as observed when injecting C particles 
into the system. Indeed, the creation of pairs of A particles on the one hand changes the noise-noise correlator which 
now reads (with h = h2A/l'^) 



(C(x, t)C{^',t')) = [h{t) - 2dt{c'^{t))] <5(x - x')5(t - t') , 
and on the other hand modifies the Langevin equation for the A particles: 

{dt ~ Da\7^)a{x, t) = -2Aa2(x, t) + 2^c(x, t) + C(x, t) + 2h{t) 



(43) 



(44) 



Using in addition that the Langevin equation for the C particles remains unchanged, it is then straightforward to 



show that we have i?^^ 



R'^ and R^a 



Re- 



The situation changes if instead of injecting pairs of A particles only single A particles are created with rate Ha- In 
that case the additional A particles do not automatically lead to the formation of additional C particles, but instead 
a newly created A particle must first diffuse through the system in order to encounter another additional A particle. 
Formally the creation of single A particles again shows up in the Langevin equation ^ for the A particles as an 



additional field term h{t), but the noise- noise correlator and, subsequently, the 6x correlator are the same as for the 
h — case: 

(C(x,i)C(x',t')) = -2ft(c"(i))<5(x-x')<5(i-i'), (45) 

{{Sx'r) - 2{co - c'^^){87rDt)-^/\ (46) 

This then readily yields the expression 



^S(i,s)-^ 



^ 2AX^i . 

1 ^ (Co - Coo) 



iTrDt)-^/^{t~s), (47) 



for the response of the C particles to the creation of single A particles. It is important to note that the expression 
p7)) does not depend (at least in leading order) on the excitation time s. This is a direct consequence of the fact that 
the rate of encounter of two newly created A particles does not depend on the time s at which the particles have been 
created. 

All the perturbations considered so far have in common that the quantity (x) = Ua + 2nc, that would be constant 
without the perturbation, has a different value in the stationary state than in the initial state. In order to assess the 
importance of this conserved quantity, we also looked at the response to a perturbation that keeps (x) unchanged. 
This can be achieved by the simultaneous creation of particles of one type and destruction of particles of the other 
type. For example, we can create pairs of A particles with rate /i,^ and at the same time remove C particles with 
rate h-^^. In principle, the removal procedure is an ill defined process, since one may end up with negative number of 
particles on a given site. However, for large equilibrium concentrations Ooo and Coo, we expect the number of particles 
on each site of the lattice to remain larger than zero following an infinitesimal excitation of the type just described. 

It is important to note that the Langevin equation (|lip for x is unaffected by this perturbation, even so additional 
field terms are entering into the Langevin equations for a and c. Doing the calculations along the lines just sketched 
for the other perturbations, one remarks that, due to the fact that the asumptotic values c^ and a^ are independent 
of the field h = h^/l'^, the response functions become time translational invariant. Thus for the response of the C 
particles to this perturbation we get 

R'^it, s) = 2^ [8TTD{t - s)]-""/^ 9{t - s) , (48) 



whereas for the response of the A particles we obtain 



i:^ra^n^+_.M-'i/2, 



i?^(i, s) = -2R^{t, s) = -4:^ [8TTD{t - s)] "I" 9{t - s) , (49) 

as il = 0. 

It has been proposed [33, [3J] that in aging systems space-time symmetries can be exploited in order to derive exact 
expressions for two-time quantities (see [6| and [7| for recent reviews of this approach in the context of magnetic and 
of reaction-diffusion systems, respectively). Thus for the autoresponse function, this theoretical approach, called the 
theory of local scale invariance, yields the following general expression: 

/,\ I+o'-Ah/z /. \ -l-a' 

R{t,s)=ros-''-^(-j (s^O ^(*-s)' (50) 

where the values of the exponents a, a' and A^^ are not fixed by the theory, whereas rg is a nonuniversal numerical 
prefactor and z is the dynamical exponent. In the reversible reaction-diffusion system considered here, we have the 
interesting situation that we can design different perturbations of the system and monitor the reaction of the system to 
these perturbations. Comparing the theoretical expression (j50p with our exact results, we observe that the responses 
(|48)) and (|49| to a perturbation that does not change the conserved quantity (x) can indeed be cast in the form (|50p . 
with a — a' = d/2 — 1 and Xr/z — d/2. On the other hand, however, perturbations that change the value of the 
conserved quantity can not be cast in the form (j50p . indicating that one must be careful when applying space-time 
symmetries to perturbations that change quantities that are otherwise conserved by the dynamics of the system. 

III. OTHER REVERSIBLE REACTION SCHEMES 
A. The A + B ^ C reaction scheme 

The treatment of the bimolecular reaction scheme A + B ^^ C closely follows that of the scheme A + A ^=^ C 
discussed in the previous Section. The main difference is the presence of two 'conservation' laws for the average 



densities |26| : 

000+^00 + 2coo = ao + bo + 2co and floo - &oo = ao - ^o , (51) 

where oq, bo, cq, respectively Ooo, &oo, Coo, are the concentrations of the A, B, C particles in the initial, respectively 
final, state. Writing down the exact Langevin equations, one notices that the Langevin equation for c(x, t) is still 
noise independent. Exploiting this property, we find for the two-point correlation function the expression 

C(t,s)-i('^') {c^-c^){A^D{t + s))-''/\ (52) 

which is of the same form as for the A + A^=^C reaction scheme, see Eq. (P^. The only difference is that the constant 
a is now given by cr = A(ao + 60 + 2co — 2coo) + jJi- We also notice that the expressions for the different response 
functions are unchanged, provided that the new expression for a is used. 

B. The A + B ^ C -\- D reaction scheme 

More interesting is the reversible reaction A + B ^^ C + D^ not discussed in [2y|, as here the Langevin equation for 
each reactant does depend on noise. It is worth mentioning that the A + B ^ C + D reaction scheme is readily found 
in experimental situations, one well-known example being ethanoic acid dissolved in water that forms ethanoate and 
hydronium ions: CH3CO2H + H2O ^^ CH3CO2 + H^O^. In order to make the following discussion more compact, 
we use the symbols Ai and A2 instead of A and B, and Ci and C2 instead of C and D. 

The exact Langevin equations for this four species reversible reaction read: 

^-Z^V^'ja, = /(a„c,) + Ca, , (53) 

^ - DvA c, = -/(a„ c) + Cc, , (54) 

where i — 1,2 and f{a,c) = (/iCiC2 — Aaia2). Here the complex variables are again related to the average particle 
densities, e.g. Ua^ — (ai) is the mean density of the particles of type Ai. The only non vanishing noise-noise 
correlators are 

(Ca,(r)Ca.(r'))=2d>-r')(/(a,c)), (55) 

(Cc,(r)Cc.(r')) = -25(r-r')(/(a,c)). (56) 

There exist three 'conserved' quantities, namely ipt = J2i=ii'^i + '^0; 4'a ^ cii — 02, and 0c = ci — C2, which obey the 
noisy diffusion equations 



DV^ <f>, = (57) 



d_ 
dt 

with j — t,a,b. The noise terms Cj are thereby just linear combinations of the CaiiCci^ Ct = X]i=i(Cai + Cci), 

Ca = Cai ~ Ca2 i and 4c = Cci ~ Cc2 • 

Let us now calculate the average density of, say, Ai particles, (ai). Since the solution of the steady state is given 
by the condition f{af^, c°°) = 0, where a°° and c°° are the average particle densities in the steady state, we can find 
all four coefficients a°° and c°° {i = 1,2) by exploiting the conservation laws {(pt) = const = $t, ((/jq) ~ const = "I>a, 
and (0c) = const = $c: 






4$f ^ 4$t 



(58) 



We will in the following discuss to some extent the case ^ — X which has the virtue that the algebra is quite easy. 
Moreover, by choosing an appropriate time rescale we can always set n = 1. The general situation A ^ /i can be 
treated along the same lines, but the algebra is rather involved. As we find for the general case at late times the 



same functional dependence of the particle density as for the case /i = A, we refrain from giving the details of this 
calculation and only quote the result at the end of this Section. 



Substituting a^ = a° 
after some algebra, 



d_ 
dt 



Odi. Ci 



DV^ 5ai 



+ &i, and <j)j — <f>j + 5(j)j (where i = 1, 2 and j — t, a, c) into Eq. (|53p yields, 



(1)2 -I- ct)2 _ $2 



4$* 



■&t>a + -Z-5(l)c + Cai 



1 



+ (l)a- 4>c) - n^t S4>a + 5(j)t 5ai 



(59) 



Note that the above equation is linear in terms of the function Sai. This is a consequence of our choice /i = A, and is 
not true in the general case. 

By introducing a new Green function 



G = 



= c-*t* 



Go(x,i), 



(60) 



we can set up a perturbation series along the same lines as discussed previously for the other reaction schemes. We 
then get 



{Sai) = -G 
+G 



l) + {&f>j)^l-{^t&l>. 



G 



$2 



^i-^i 



^t~^a 



4$t 



-^ 5/)c + Cai 



(61) 



where the dots refer to higher order terms that are omitted from now on. Luckily, only few of the many terms in this 
equation are different from zero. Indeed, since {&t)a{r) S(f)a{v')) = —{&j)a{^)^a{^')) and {&t)a{r) 5cf>c{r')) = 0, we get 
{Scj)t{r) §cj)t{r')) — and {S(f>t{r) 6ij)j{r')) — where j = a, c. This follows directly from the equations (fSS]) and (fSG]) . It 
is also clear that the noise term contribution (^ai in Eq. (|6ip is exponentially suppressed and hence can be discarded. 
With this we end up with the following expression: 



(<5ai) 



-G 



\{&f>i-m 



= 1^ y" dti y" d-^xi G2(x - xi, t - ii)(/(a„ cO(ii)) . 



(62) 



from which wc obtain that the approach to the equilibrium density of Ai particles is again governed by a power law, 
namely 



(Sai) = — [ai(oo) - ai(0)](87rDt) 



-d/2 



(63) 



Following the same stepe as in Section II, we obtain for the autocorrelation function the expression 

($t - $a)' - *? 



C{t,s) 



4$2 



■-{&j)a{t)&j)a{s)) 



<^i 



2<D2 



[ai(oo)-ai(0)][47r7:>(i + s)] 



-d/2 



(64) 



with the same dependence on s and t as for the A + A ^ G reaction scheme. Response functions are also calculated 
as previously, yielding again the same functional dependences. For example, for the response of the Ai particles to a 
perturbation that conserves <&< we get 



Ri: (*' ^) = :5- [«i (°o) - «i (0)] [8^-0(i - s)] -"!'' 9(t - s) . 



^t 



(65) 



Let us finish this Section by very briefly discussing the general case /i 7^ A. In principle, one follows exactly the 
same steps as before, but the algebra is much more involved. Keeping only relevant terms we get instead of Eq. ((59|) 
the expression 



d 
dt' 


J 


Sai = 


- (Tdai 2 «j)a+ ^ 



(A — ij)&if 



(66) 
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with a = -^//zA^j + (A — /x)(A$^ — /i$c)- This then yields the following asymptotic approach toward the equilibrium 
density of the A\ particles: 



((fai) = g(Ai, A, $„ $,)[ai(oo) - ai(0)](87ri?t)-'^/2 , (67) 



where we defined 



5(/i, A, $,, $,) = ^^ + ^ [(3a - A$,)2 - (^$,)2] . (68) 

It is again to be noted that this leads to exactly the same functional forms for the autocorrelation and autoresponse 
functions. 

IV. CONCLUSION 

In this work we extended the study of aging phenomena in reaction-diffusion systems toward reversible reaction 
schemes. Our starting point thereby was the observation [1^ that in reversible reaction-diffusion systems slow dy- 
namics, i.e., dynamics characterized by power law relaxation, is generic. This is in stark contrast to the case of 
irreversible reactions as here the approach to stationarity usually happens exponentially fast, with the exception of 
nonequilibrium critical points which are governed by power laws. 

In our study we focused on simple models with reversible reactions whose behavior in the asymptotic or scaling 
regime can be computed exactly j26|. We observed for all studied models a simple scaling behavior of two-time 
response and correlation functions. Interestingly, and in agreement with irreversible reaction-diffusion systems at 
their critical point, the two scaling exponents a and b, see Equations ((2]) and ([T]), were found to be different. The 
multi-species models we studied have the virtue that we can define different responses, depending on the particle type 
for which an additional creation process is considered and on the particle type whose density we are monitoring after 
the perturbation of the system. We also studied a response where the perturbation conserves the total number of 
particles in the system. For this perturbation we found that the exact result agrees with the expression one obtains 
from the theory of local-scale invariance [y, [33f| that exploits the existence of space-time symmetries in aging systems. 
However, if the system is perturbed in such a way that the total number of particles, a quantity that is constant in 
the unperturbed system, is not conserved, than the response can not longer be cast in the theoretically predicted form 
(j50p . This result indicates that some care has to be taken if one wants to apply space-time symmetries to cases where 
response functions result from perturbations that change some otherwise conserved quantities. 

It is very appealin g th at the theoretically predicted power law approach to stationarity is readily observed in 
experimental systems [30, |3l|, |33 • It should therefore be possible to measure in these systems two-time quantities in 
the scaling regime in order to verify the scenario of simple aging that follows from our study. 

The models we consider in this work are to some extend artificial, as we allow multiple occupancy of a given site and 
only consider on-site reactions. It is an important question whether the scaling picture emerging from the study of 
these simple models also holds in more realistic cases. Of special interest in this context is the restriction to single-site 
reactions, as this disagrees with the actual experiments where longer-range reactions prevail. Indeed, it has been 
stressed in the literature 26] that the models studied in this paper yield as a stationary state a chemical equilibrium 
state. Voituriez et al. [35], however, pointed out that distance-dependent reversible reaction rates no longer yield 
asymptotically a chemical equilibrium state, but that the stationary state is then a nonequilibrium state. This raises 
the interesting prospect that one could be able to study the similarities and differences in the aging behavior of 
systems relaxing toward equilibrium and nonequilibrium stationary states by changing the range of the reactions. We 
plan to study this important aspect in our future work. 

APPENDIX A: RELATION BETWEEN THE ASYMPTOTIC PARTICLE CONCENTRATIONS 

In the discussion of the aging processes for the A + A ^ C reaction we extensively used the relation pH)) between 
the equilibrium densities of the A and C particles. In this Appendix we briefly show that in zero dimension this 
condition may not be realized by all initial conditions. 

Starting from the probability P(m, n, t) for having at time t m resp. n particles of type A resp. C at some lattice 
point, we can introduce the generating function (see, for example, |lll.l3a|) 

^{d,d,t) = Y,a"'c"P{m,n,t). (Al) 
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It then follows from the Master Equation (|4]) that this generating function is solution of the partial differential equation 

It is more convenient to work with the shifted variables a = a — 1 and c = c — 1 which yields for the equation for the 
generating function the expression 

|*.(a" + 2S-o-)(4-A^)*, (A3) 

The stationary solutions of this equation are given by 

/oo 
dzGn{c,a~z)f{z), (A4) 

-CxD 

where Go is just the Green function for the one-dimensional diffusion equation, 

y , s -1/2 / 2 \ 

Go(x,y) = Ut:-x] exp I -^ J , (A5) 

and f{z) is any function satisfying the condition /(O) = 1 with only positive coefficients in the Taylor expansion. It 
directly follows that there exist infinitely many solutions of the equation (jASp . such as 

*o(a,c) = l + -c + yla+-a^ (A6) 

(with an arbitrary constant A) or 

'A 



*p(a, c) = exp ( -c + a 1 , (A7) 

which corresponds to the Poisson stationary distribution. It is now easy to see that the condition iiCoo = Aa^ is 
indeed realized for the Poisson distribution (IA7D. but not for the distribution (IA6I). 



APPENDIX B: PERTURBATION SERIES 

In their article [2^ Rey and Cardy stress that one of the central points of their analysis is the fact that the equation 
for the variable c comes without explicit noise, thus implying that their approach only works in cases where one of 
the equations of motion does not have an explicit noise dependence. In fact, this requirement is not really needed, 
as we demonstrate in this Appendix by deriving the average number of A particles directly from the 'noisy' equation 

Indeed, starting from this equation we have 

(dt ~ DaV^ + ^a)Sa = ^fi5x + C - A&i^ . (Bl) 



AG [(fa2] , (B2) 



The formal solution of the nonlinear equation (JB1[) is given by 

'1 



Sa^G 



2^'^ + C 



where the Green function G is the inverse of the operator dt — DgS/^ + ^o'- With this one can derive a perturbation 
series for, say, the average density {5a) following the same strategy as in [2a|. One thereby exploits the fact that the 
noise term contribution G[Q in this series is exponentially suppressed and hence can be discarded. To convince oneself 
that it is indeed so, one should consider terms such as (G[C]^) and (G[(5\;]G[C]), which are exponentially small for late 
times t. Taking this into account, one gets 

{5a) = X (lA "' a^XiSx') = -^-^ie) , (B3) 
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which, of course, corresponds to the result 

{Sa) = -2(Sc) . (B4) 

obtained in [2^| as the condition {6x) = holds. 
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